A light-cone gauge for black-hole perturbation theory 
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The geometrical meaning of the Eddington-Finkelstein coordinates of Schwarzschild spacetime 
is well understood: (i) the advanced-time coordinate v is constant on incoming light cones that 
converge toward r = 0, (ii) the angles 8 and <f> are constant on the null generators of each light cone, 
(iii) the radial coordinate r is an afnne-parameter distance along each generator, and (iv) r is an 
areal radius, in the sense that inr 2 is the area of each two-surface (v,r) — constant. The light-cone 
gauge of black-hole perturbation theory, which is formulated in this paper, places conditions on 
a perturbation of the Schwarzschild metric that ensure that properties (i)-(iii) of the coordinates 
are preserved in the perturbed spacetime. Property (iv) is lost in general, but it is retained in 
exceptional situations that are identified in this paper. Unlike other popular choices of gauge, the 
light-cone gauge produces a perturbed metric that is expressed in a meaningful coordinate system; 
this is a considerable asset that greatly facilitates the task of extracting physical consequences. We 
illustrate the use of the light-cone gauge by calculating the metric of a black hole immersed in a 
uniform magnetic field. We construct a three-parameter family of solutions to the perturbative 
Einstein-Maxwell equations and argue that it is applicable to a broader range of physical situations 
than the exact, two-parameter Schwarzschild-Melvin family. 
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I. INTRODUCTION 

The theory of perturbations of the Schwarzschild 
spacetime is a well-developed one 0], and it may seem 
surprising that authors are still writing on this venerable 
topic almost fifty years after its inception in the work 
of Regge and Wheeler (see also Refs. 0,01). This is 
(at least partially) explained by the fact that the field 
has witnessed a resurgence of sorts in the last several 
years, motivated by new applications that include the 
gravitational self- force problem 0, 0j the "close- 
limit" collision of two black holes 0], and the study of 
the dynamics of black holes placed in tidal environments 
M Ell IH Elll The theory has been presented in 
various sophisticated packages 

m M El El MM E3, 

and it has reached what is likely to be its definitive form. 

We wish to make an additional contribution to this 
body of literature by formulating a useful, attractive, and 
simple gauge condition for black-hole perturbation the- 
ory. We believe that this gauge, which we call the light- 
cone gauge, is preferable in many ways to most popular 
gauges, including the oft-used Regge- Wheeler gauge. We 
believe that the use of the light-cone gauge will be a great 
benefit to any researcher faced with the task of comput- 
ing and interpreting a perturbation of the Schwarzschild 
spacetime. 

The idea is simple. The difficulties of the Schwarzschild 
coordinates (t, r, 9, 4>) across the black-hole horizon are 
well documented, and it is well known that the transfor- 
mation v = t + r + 2M\n(r/2Al — 1) brings the metric 
to a form that is well-behaved on the event horizon. The 
Eddington-Finkelstein coordinates (v, t, 9, cf>) have a clear 
geometrical meaning. The null coordinate v (called ad- 
vanced time) is constant on incoming light cones that 
converge toward r — 0, the angles 9 and <f> are constant 
on the null generators of each light cone v — constant, 
and r is an afnne-parameter distance along each genera- 



tor. In addition, r doubles as an areal radius, in the sense 
that 4irr 2 is the area of each two-sphere (v, r) — constant. 

The light-cone gauge places conditions on the metric 
perturbation that ensure that the geometrical meaning of 
the coordinates is preserved in the perturbed spacetime. 
The advanced-time coordinate v therefore continues to 
label incoming light cones that converge toward r = 0, 
the angles 9 and <j> continue to label the generators of each 
light cone, and r continues to be an afnne-parameter dis- 
tance along each generator. One geometrical aspect of 
the coordinates that must generally be given up is the 
role of r as an areal radius; we shall show, however, 
that this property also can be preserved in special cir- 
cumstances. The light-cone gauge therefore produces a 
perturbed metric that is expressed in a meaningful coor- 
dinate system. This is a considerable asset that greatly 
facilitates the task of extracting the physical properties 
of the spacetime. 

The light-cone gauge is developed in Sec. II of this 
paper. The gauge conditions are introduced in Sec. II 
B, after we present in Sec. II A a brief review of the 
Schwarzschild metric in Eddington-Finkelstein coordi- 
nates. In Sec. II C and D we decompose the metric per- 
turbation in spherical harmonics and explore the space of 
gauge transformations that keep the perturbation within 
the light-cone gauge. This remaining gauge freedom is 
convenient, as it can be exploited to simplify the form 
of the perturbed metric to the fullest extent possible. 
In Sec. II E we determine the conditions under which 
r retains its interpretation as an areal radius. The an- 
swer turns out to be simple: This is possible whenever 
the component T vv = T rr of the perturbing energy- 
momentum tensor vanishes. In Sec. II F we summarize 
our construction and discuss its merits; in particular we 
compare our light-cone gauge to the very widely used, 
but far less compelling, Regge- Wheeler gauge. 

In Sec. Ill we present an illustrating application of the 
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light-cone gauge for black-hole perturbation theory: We 
examine a black hole immersed in a uniform magnetic 
field, and calculate its metric accurately through second 
order in the strength of the magnetic field. The physical 
situation is described in Sec. Ill A. The magnetic field 
and its energy-momentum tensor are computed in Sec. Ill 
B. In Sec. Ill C and D we integrate the equations of 
black-hole perturbation theory for this situation. The 
solution is presented in Sec. Ill E, and in Sec. Ill F we 
examine the structure of the perturbed horizon. Finally, 
in Sec. Ill G we compare our perturbative solution to the 
exact Schwarzschild-Melvin solution [22l I23I |24| , which 
also describes a magnetized black hole. We conclude that 
the perturbative solution is applicable to a broader range 
of physical situations. 

In the Appendix we provide a complete listing of the 
linearized field equations in the light-cone gauge. 

Our developments in this paper rely heavily on the re- 
cent work of Martel & Poisson [2l| , which presents a co- 
variant and gauge-invariant formalism for black-hole per- 
turbation theory. They also incorporate key ideas from 
a companion paper ;25i] devoted to the construction of 
light-cone coordinates centered on a geodesic world line 
of an arbitrary curved spacetime. 

We point out that the light-cone gauge constructed 
here is adapted specifically to incoming light cones v = 
constant that converge toward r — 0. It would be exceed- 
ingly straightforward to adapt the construction to outgo- 
ing light cones u = constant that expand toward r — 00. 
[In Schwarzschild spacetime, the retarded time coordi- 
nate u is defined by u = t-r-2Mln(r/2M- 1).] While 
the incoming light-cone gauge is well suited to study the 
properties of the perturbed event horizon, the outgoing 
light-cone gauge would be well suited to study the grav- 
itational radiation escaping toward future null infinity. 
We suggest this adaptation as an exercise for the reader. 

Throughout the paper we work in geometrized units 
(c = G = 1) and adhere to the conventions of Misner, 
Thorne, and Wheeler 12611. 



II. LIGHT-CONE GAUGE: DEFINITION AND 
PROPERTIES 

A. Schwarzschild metric in light-cone coordinates 

The transformation v = t + r + 2Mln(r/2M — 1) 
brings the Schwarzschild metric from its usual form to 
the Eddington-Finkelstein form 



ds z = -f dv z + 2 dvdr + r 2 dti 2 



where 



/:=1-™ 
r 



(2.1) 



(2.2) 



and 



The parameter M is the black-hole mass and 9 A — 
(9 2 ,9 3 ) — (9,4>) are angles that span the two-spheres 
(v, r) — constant. 

The coordinates (v, r, 9, (j>) are well-behaved across the 
event horizon, and they possess a clear geometrical mean- 
ing. We note first that the vector 



l a :=-V a v = (-1,0,0,0) 



(2.4) 



is null, which implies that each surface v = constant is a 
null hypersurface of the Schwarzschild spacetime; these 
are in fact incoming light cones that converge toward 
the black-hole singularity. The fact that l a is a gradient 
implies that l a is everywhere tangent to a congruence of 
null geodesies; these are affinely parameterized, and they 
are the generators of each light cone v = constant. Using 
the metric of Eq. 1|2.1[) to raise indices, we find that 



r = (0,-1,0,0). 



(2.5) 



This relation implies that 9 A = constant on the gener- 
ators, so that the angles 9 A can be used as generator 
labels. Furthermore, Eq. (|2.5|l reveals that the affine pa- 
rameter on each generator is — r. The geometrical mean- 
ing of the coordinates is therefore the following: The null 
coordinate v (called advanced time) is constant on incom- 
ing light cones that converge toward r = 0, 9 A labels the 
generators of each light cone, and r is an affine-parametcr 
distance along each generator. The radial coordinate r 
also doubles as an areal radius, meaning that 4irr 2 is the 
area of each two-sphere (v, r) = constant. 



B. Perturbed metric in light-cone coordinates 

We introduce a perturbation p a p of the Schwarzschild 
metric, defined by the statement 



perturbed 
9a/3 =g a p+Paf3, 



(2.6) 



dQ 2 := n AB d9 A d9 B := d9 2 + sm z 9 d<f>. (2.3) 



where <?^ rturbcd is the metric of the perturbed spacetime, 
while g a p represents the Schwarzschild solution, which 
we express in the light-cone coordinates of Eq. (|2.1|) . We 
wish to place conditions on the metric perturbation that 
ensure that the meaning of the light- cone coordinates will 
be preserved in the perturbed spacetime. Specifically, we 
demand that in the perturbed spacetime, v continues to 
be constant on incoming light cones that converge toward 
r = 0, 9 A continue to be constant on the null genera- 
tors of each light cone, and r continues to be an affinc- 
parameter distance along each generator. In exceptional 
circumstances that will be identified in Sec. II E below, 
r also reprises its role as an areal radius, but in gen- 
eral this property will not be preserved in the perturbed 
spacetime. 

The geometrical meaning of the coordinates will be 
preserved if Eqs. I|2.4|l and l|2.5|l continue to hold in the 
perturbed spacetime. We now have l a = (g a p + Pap)! 13 , 
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and we infer that the perturbation must satisfy the gauge 
conditions 



Papl = 



PrA 



0. 



(2.7) 



There are four conditions, which we refer to as the light- 
cone gauge conditions. A metric perturbation p a p that 
satisfies Eqs. (|2.7|) will be said to be in a light-cone gauge. 
As we shall see below, Eqs. (|2.7|l do not completely fix the 
gauge, and the remaining gauge freedom can be exploited 
to simplify the form of the perturbed metric. 

The gauge conditions leave p vv , p v A, and pab as non- 
vanishing components of the metric perturbation. The 
trace of the perturbation is 



p := 



g aP Pa0 



r-*n AB p AB , 



(2.8) 



where il AB is the matrix inverse of fi^, defined by 
Eq. The determinant of the perturbed metric is 

given by 



v /_ g perturbed _ ^ sm ^l _|_ _ p 



(2.9) 



and r will retain its role as areal radius whenever the 
metric perturbation has a vanishing trace. In Sec. II E 
we will determine under what conditions this happens. 



C. Even- parity sector 

The even-parity sector refers to those components of 
the metric perturbation that can be expanded in terms 
of even-parity spherical harmonics Y lm , Y A m , VIabY 1 ™ 1 , 
and F^g. (Throughout the paper we use the notation 
of Martel & Poisson |2l|.) The scalar harmonics Y lm 
are the usual spherical-harmonic functions, the vecto- 
rial harmonics are defined by Y A m := D aY 1 ™ (where 
Da is the covariant derivative operator compatible with 
QjIb), and the tensorial harmonics are defined by Y A r g := 
[DaDb + + 1)£Iab] Y lm ; these are tracefree by virtue 
of the eigenvalue equation for the spherical harmonics: 

ClAByhn = \£l AB D A D B + 1(1 + l)]Y lm = 0. 

The even-parity sector is 



Pab 
PaB 
PAB 



^C(^)F ;m (^), (2.10) 

Irn 

J2j a m (x a )Y B m (6 A ), (2.11) 

Irn 

r 2 [K lm {x a )n A BY lm {e A ) 



I in 



+ G lm (x a )Y A n B(9 A ) , (2.12) 



where x a ~ (x ,^ 1 ) = (v,r). The sums over the integer 
I begin at I = in the case of Eq. (|2.1U|I and the K lm 
term in Eq. (|2.12|) . at I = 1 in the case of Eq. I|2.11|l . and 
at I = 2 in the case of the G lm term in Eq. $FF^i . The 



sums over the integer m go from — I to +1. The light-cone 
gauge conditions are 



Ira 



,trn _ ..... = „, 

' "VT TT Jy 



Irn 



0. 



(2.13) 



The components h l ™, j l v m , K lm , and G lm are nonzero in 
the light-cone gauge. 

Even-parity gauge transformations are generated by a 
dual vector field S Q = (S a , S^) that can be expanded as 



Irn 

za = E^ m (^)^ 



Imi nA 



A ). 



(2.14) 
(2.15) 



Irn 



According to Eqs. (4.6)-(4.9) of Martel & Poisson |21j . 
such a transformation will preserve the conditions of 
Eq. (|2.13|l provided that £„, £ r , and £ satisfy the equa- 
tions 



dr dv 
dr ' 



2M 



dr 



&■ - -e 



(We henceforth omit the spherical-harmonic indices for 
brevity. Our considerations momentarily exclude the 
special cases I = and 1 = 1, which will be handled 
separately below.) This means that a gauge transforma- 
tion generated by 

& = -a(v)r-fa(v) + b(v), (2.16) 
Hr = a(v), (2.17) 
£ = a(v)r + c(v)r 2 (2.18) 

will keep a perturbation within the light-cone gauge. The 
remaining gauge freedom is therefore characterized by 
three arbitrary functions a(v), b(v), c(v), and the over- 
dot in Eq. I|2.16|l indicates differentiation with respect to 
v . The gauge transformation changes the nonvanishing 
components of the perturbation field according to 



J 3M\. 
= h vv + 2a r + 2 1 a 



Jv * Jv 

K K' = K 



2M, 

jv + fa 



2a 



b- 
1(1 



r J 



cr 2 , 

+ 1) 2 
a b 

r r 



2b 



+ l(l + l)c, 



G^G' = G a - 2c. 

r 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



The lower multipoles I = and I = 1 must be consid- 
ered separately. For I = the spherical harmonics Ya 
and Yab are identically zero, and the only relevant per- 
turbation fields are h a b and K; j a and G are not defined. 
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A gauge transformation is then generated by S a = £, a Y 00 , 
"E, A = 0, and £ is not denned. It is easy to check that the 
light-cone gauge will be preserved with a £ a that is still 
given by Eqs. I|2.16[l and l|2.17[) . In this case the remain- 
ing gauge freedom is characterized by two arbitrary func- 
tions, a(v) and b(v). The corresponding change in h vv is 
still given by Eq. i|2.19[l . while K' is obtained by setting 
I = in Eq. (|2.21|) . For I = 1 the tensorial harmonics 
Yab are identically zero, and only G is not defined. The 
gauge transformation of Eqs. (|2.16l) - (|2.18|) is still seen to 
preserve the light-cone gauge, and the changes in the per- 
turbation fields are still described by Eqs. H2.19fl - 12.21fl . 
in which we must set 1(1 + 1) = 2; Eq. (|2~2l|) is irrelevant 
when 1 = 1. 

It is easy to verify that the dual vector field of 
Eqs. (|2.14l) - (|2.18|) generates the (small) coordinate trans- 
formation 



v' = v + a(v,9 A ), 
da 



nlA 



1 



= e J 



dv 



r + b{v,9 A ), 



AB 



d 



+ c{v,9 A ) 



(2.23) 
(2.24) 

,(2.25) 



where a(v, 9 J 



J2 lm a lm (v)Y lm (6 A ), with similar 
This transfor- 
mation leaves the conditions of Eq. I|2.7|l intact. 



equations defining b(v,9 ) and c(v,9 



D. Odd-parity sector 



The odd-parity sector refers to those components of 
the metric perturbation that can be expanded in terms 

and X%. The 

-efD B Y lm , 

where eab is the Levi-Civita tensor on the unit two- 
sphere (with independent component ee^ = sin 9), and 

yBC = 



of odd-parity spherical harmonics X l Y' 
vectorial harmonics are defined by X A :- 



where e A := VI^^eac- The tensorial harmonics are 
:= -\{e£D B +sgD A )D c Y lr ' 



AB 



these are tracefree 



by virtue of the antisymmetry of the Levi-Civita tensor: 



CtAB ylm 

11 A AB — 



e AB D A D B Y lm = 0. 



The odd-parity sector is 
Pab = 0, 



PaB = Y,^{x a )X l B{0 A ), 



In i 



PAB 



^h l 2 m (x a )X%(9 A ). 



I III 



(2.26) 
(2.27) 

(2.28) 



The sums over the integer I begin at I = 1 in the case of 
Eq. (1223, and at I = 2 in the case of Eq. The 
light-cone gauge conditions are 



Am 



0. 



(2.29) 



The components h l ™ and h ™ are nonzero in the light- 
cone gauge. 



Odd-parity gauge transformations are generated by a 
dual vector field S a = (H ,H^) that can be expanded as 



S a = 0, 



-.4 



I in 



C tm (x a )X l A m (9 A ). 



(2.30) 
(2.31) 



According to Eqs. (5.5) of Martel & Poisson [2l(, such a 
transformation will preserve the conditions of Eq. (|2.29(l 
provided that £ satisfies d^/dr — 2£/r = 0. (We hence- 
forth omit the spherical-harmonic indices for brevity. 
Our considerations momentarily exclude the special case 
1 = 1, which will be handled separately below.) This 
means that a gauge transformation generated by 



£ = a(v)r 2 



(2.32) 



will keep a perturbation within the light-cone gauge. The 
remaining gauge freedom is therefore characterized by a 
single arbitrary function a(v). The gauge transformation 
changes the nonvanishing components of the perturba- 
tion field according to 



h v 
h 2 



h' v = h v — dr 2 , 
ti 2 = h 2 - 2ar 2 , 



(2.33) 
(2.34) 



where an overdot indicates differentiation with respect to 
v. 

The situation is the same for the special case / = 1, 
except that X A b is then identically zero and hi is not 
defined. The gauge transformation of Eq. I|2.32|l is still 
seen to preserve the light-cone gauge, and it still changes 
h v according to Eq. I|2.33|l : Eq. I|2.34|l is then irrelevant. 

It is easy to verify that the dual vector field of 
Eqs. (|2.30(l - (|2.32|l generates the (small) coordinate trans- 
formation 



AB 



d_ 

89^ 



a(v,9 A ), 



(2.35) 



where a(v, 9 A ) 



(v)Y im (9 A ) and 



-AB 



Ecd- This transformation leaves the conditions 
of Eq. lt2~7l) intact. 

E. When is r an areal radius? 

According to Eqs. (|2.8|l and (|2.9() . r keeps its interpre- 
tation as an areal radius whenever £I ab pab = 0. And 
according to Eq. fHfy . this happens when K lm (v, r) = 
for all values of I and m. In this subsection we determine 
under what circumstances it is possible to impose this 
condition. 

The light-cone gauge produces a very convenient de- 
coupling of the equation that governs the behavior of K lm 
from the equations that determine the remaining pertur- 
bation fields. According to the field equations listed in 
the Appendix, we have 



trn 



dr 2 



K 



lm 



2d_ 

r dr 



K 



lm 



(2.36) 
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where, for example, Q l ™ := 8tt J T rr r im dO, with dfi 
si 1 1 dOdcj), are the spherical-harmonic projections of the 
rr component of the energy-momentum tensor. When 
T vv = T rr = 0, Eq. QHty reveals that K lm must be 
of the form p lm (v) + q lm (v)/r, where p lm and q lm are 
arbitrary functions of v. But it is possible to exploit 
the remaining gauge freedom contained in Eqs. I|2.16[l - 
(|!Tr%|) to set K" m = 0. As can be seen from Eq. QSD , 
this condition constrains the functions b lm (v) and c lm (v), 
which must now be related to a lm {v). The remaining 
gauge freedom is therefore restricted to transformations 
characterized by a single arbitrary function, a lm (v). Our 
conclusion is that K lm can be set equal to zero whenever 
rpvv _ rp rr _ anc j that this operation still does not 
fully exhaust the gauge freedom. 

We have established the following theorem: When the 
energy-momentum tensor responsible for the metric per- 
turbation is such that 

T aP l a f = 0, (2.37) 

the light-cone gauge can be refined to include the trace- 
free condition 

P ■= g aP p af 3 = (2.38) 

in addition to the four conditions of Eq. I|2.7|) . In such 
circumstances, \/— 3 pcrtulbcd = r 2 sin (9 and r retains its 
interpretation as an areal radius. In these circumstances 
the light-cone gauge becomes the "incoming radiation 
gauge" of Chrzanowski H |H El HI • 

The theorem relating the tracefree condition of 
Eq. (|2.38|) to the vanishing of T a pl a lP is a new result. 
The theorem was established independently by Price, 
Shankar, and Whiting 0, m work that has not yet 
been published, except for a statement of the result made 
in Sec. 4.1 of Ref. [33j. Remarkably, these authors were 
able to extend the theorem from Schwarzschild spacetime 
to all Petrov type-II spacetimes. 

F. Discussion; Comparison with the 
Regge- Wheeler gauge 

The light-cone gauge possesses two main virtues. The 
first is that it involves simple algebraic conditions on the 
metric perturbation; these were stated in covariant form 
in Eq. I|2.7|l . p a pl^ — 0, and they were stated in expanded 
form in Eqs. lf!TT3j) and h vr = h l"r = J 1 " 1 = K™ = 

0. The second is that the gauge conditions preserve the 
geometrical meaning of the original coordinate system 
(v, r, 6 ); as was shown in Sec. II B, the advanced-time 
coordinate v continues to label incoming light cones that 
converge toward r = 0, the angles 6 A continue to label 
the generators of each light cone, and the radial coordi- 
nate r continues to be an affine-parameter distance along 
each generator. The task of extracting the physical prop- 
erties of a perturbed spacetime will be greatly facilitated 
by the use of such meaningful coordinates. 



Most of the literature on black-hole perturbation the- 
ory employs an alternative gauge known as the "Regge- 
Wheeler gauge" . The gauge conditions in this case are 
3v n = fr m = G lm = hip = 0. The Regge- Wheeler gauge 
also has the advantage of involving simple algebraic con- 
ditions on the metric perturbation. But unlike the light- 
cone gauge, the Regge- Wheeler gauge produces a coor- 
dinate system that does not possess a clear geometrical 
meaning; this is a disadvantage. And indeed, the coordi- 
nates can sometimes be pathological. For example, the 
Regge- Wheeler gauge produces metric components that 
do not display asymptotically-flat behavior near future 
null infinity, even when the source of the perturbation is 
spatially bounded Q. This problem is associated with 
the fact that by imposing G lrn = h 1 ™ = 0, the Regge- 
Wheeler gauge is actually setting to zero the transverse- 
tracefree part of the metric perturbation, thereby ef- 
fectively "gauging away" its gravitational- wave content. 
(The gravitational-wave modes are still present in the 
metric perturbation, but in the Regge- Wheeler gauge 
they are encoded in unnatural places.) The end result 
is a meaningless coordinate system, a metric perturba- 
tion that fails to be asymptotically flat, and a spacetime 
that does not easily reveal its radiation content. These 
problems are not present in the light-cone gauge. 

Another approach that has been followed in the litera- 
ture on black-hole perturbation theory is to avoid fixing 
the gauge, and to work instead with a gauge-invariant 
formalism HI H El E E M IHl . Such an approach 
can be very useful, especially when an application calls 
for a switch from one gauge to another. We would argue, 
however, that a good choice of gauge can be even more 
useful in concrete applications. After all, most relativists 
would begin an investigation of the Schwarzschild space- 
time by making a specific choice of coordinate system; 
few relativists would insist on staying coordinate free. 
And most relativists would agree that the Eddington- 
Finkelstein system (v,r,9 A ) is more convenient to work 
with than the Schwarzschild coordinates (t,r, 9 A ) when 
one is concerned with the event horizon; these relativists 
would say that the Eddington-Finkelstein coordinates are 
good coordinates. These attitudes need not change when 
one goes slightly away from the Schwarzschild spacetime, 
and the light-cone gauge provides a good coordinate sys- 
tem to investigate the perturbed spacetimes. 



III. BLACK HOLE IN A MAGNETIC FIELD 

A. Physical situation 

To illustrate the use of the light-cone gauge in black- 
hole perturbation theory, we work through a model prob- 
lem involving a black hole immersed in a uniform mag- 
netic field. We have in mind a situation in which a large 
mechanical structure, such as a giant solenoid, is set up 
in outer space and made to produce a uniform magnetic 
field of strength B. The structure has a mass M' and its 
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linear extension is of the order of the length scale 1 a; the 
magnetic field is imagined to be uniform over a region 
of this size. A black hole of mass M, initially isolated, 
is then brought to the structure and inserted within the 
magnetic field. This process is quasi-static and reversible, 
and the black hole's surface area stays constant during 
the immersion. We wish to study how the black hole dis- 
torts the magnetic field within the structure, and how the 
magnetic field distorts the geometry of the black hole. 

We suppose that the perturbation created by the mag- 
netic field is small and that its effects can be adequately 
calculated with the equations of black-hole perturbation 
theory. We shall see below that the criterion for this is 
r 2 B 2 <C 1, where r is the distance from the black hole. 
If we restrict our attention to the interior of the mechan- 
ical structure and impose the inequality r < a, then the 
perturbative criterion is 



a 2 B 2 < 1. 



(3.1) 



In addition to Eq. (|3.1[1 we assume that the structure is 
situated in the black hole's weak-field region, so that 



M 

— <1. 

o 



(3.2) 



While a 2 B 2 and M/a must both be small, their relative 
sizes are not constrained. We may imagine that M/a 
is either much smaller than, comparable to, or much 
larger than a 2 B 2 ; black-hole perturbation theory can 
handle all these situations. Below we will be particularly 
(but not exclusively) interested in the first possibility, 
M/a < a 2 B 2 or M/a 3 < B 2 . In this situation there 
exists an asymptotic region (described by M < r < a) 
in which the gravitational effects of the magnetic field, 
though small, are larger than those of the black hole. 

Another aspect of our model problem is the tidal grav- 
ity exerted by the mechanical structure. Because the 
structure has a mass M' and is situated at a distance a 
from the black hole, the tidal field (or Weyl curvature) it 
produces near the black hole is £ ~ M'/a 3 . This quan- 
tity £ , which will be formally introduced below, is an 
additional parameter that characterizes the physical sit- 
uation. Below we will imagine that £ is of the same order 
of magnitude as B 2 , so that M'/a 3 ~ B 2 . Our results, 
however, will not be tied by this assumption; they will 
be just as valid when £ is much smaller than (or indeed 
much larger than) B 2 . 

The perturbed black-hole solution that we construct 
below is in fact a three-parameter family of solutions; 
each solution is characterized by the black-hole mass M, 
the magnetic field strength B, and the tidal gravity £ . 
The solution is obtained perturbatively through order 
(B 2 ,£). 



B. Magnetic field 

We first calculate the electromagnetic field F a p that 
surrounds the black hole. Because we seek to determine 
the perturbed metric accurately through order B 2 , it is 
sufficient to calculate F a p to order B. And because the 
metric corrections of order B 2 do not enter this calcu- 
lation, we may let the spacetime have an unperturbed 
Schwarzschild metric. 

To find the electromagnetic field we rely on Wald's 
observation (.3 11 ] that in a vacuum spacetime, any Killing 
vector can be identified with the vector potential of a test 
electromagnetic field. The fact that the vector satisfies 
Killing's equation ensures that the resulting F a p satisfies 
the sourcefree Maxwell equations. To produce a magnetic 
field that is asymptotically uniform when r> M, we set 



.4° 



(3.3) 



where <fi a :— (0,0,0,1) is the rotational Killing vector 
of the unperturbed Schwarzschild spacetime; we use the 
ordering (v, r, 0, (f) of the unperturbed light-cone coordi- 
nates. 

The vector potential gives rise to an electromagnetic 
field tensor F a p — VqA^ — \7pA a . To display its com- 
ponents it is useful to decompose it in an orthonor- 
mal tetrad e" that is oriented along the "Cartesian 
directions" associated with the "spherical coordinates" 
(r, 0, </>); here the superscript a is the usual vectorial in- 
dex, and the subscript [i is a frame index that identifies 
each member of the tetrad. We thus introduce the tetrad 

e£ = (r 1/2 , 0,0,0), (3.4) 
e" = (/ _1 ' 2 sin0cos<?!>, J 1 / 2 sin cos 0, r' 1 cos cos 0, 

— r _1 sin</>/sin0), (3.5) 
e 2 = {f 1 ^ 2 sin sin 0, f 1 ^ 2 sin sin 0, r _1 cos sin 0, 

r _1 cos 0/ sin 0), (3.6) 
e% = (f- 1/2 cos 9, f 1/2 cos 0, -r" 1 sin 0,0). (3.7) 



1 The constant a is not to be confused with the functions a lm (v) 
introduced in Sec. II C. 



We may think of e" as pointing in the "x-direction," of ef 
as pointing in the "y-direction," and of e" as pointing in 
the "^-direction." In this tetrad, the nonvanishing frame 
components of the electromagnetic field tensor are 

B\ := F-23 := -F^ef ef 

= v7)sin0cos0cos<£, (3.8) 

Bi := F 3 i := F a pe 3 e 1 

sin cos sin 0, (3.9) 
B 3 := F\2 '■= F a peie 2 

= B[v7 +(1-V7) cos 2 0], (3.10) 

where / := 1 — 2M/r. The field is purely magnetic. 
Asymptotically, when r > M, B x ~ 0, B 2 ~ 0, B 3 ~ B; 
the magnetic field is uniform and aligned with the z- 
direction. Closer to the black hole the magnetic field is 



7 



B\ 



distorted; at r = 2M we have \B\ 2 

B 2 cos 2 9, which indicates that the field is strongest at 
the poles. 

The electromagnetic field produces an energy- 
momentum tensor given by 



_L / p°apP 

47T * 



1 



7.9 



7<5 



(3.11) 



(3.12) 



Its nonvanishing components are 

p>2 
47T 

T vr = £-i[ r -2M-2(r-M) cos 2 0], (3.13) 
B 2 r — 2M 

T rr = s — [r-2M-2(r-M)cos 2 6>],(3.14) 



8tt r 2 

T 1 "* = sin 6 cos 0, 

4-7T r 



B 2 r — 2M 
47r r 2 

87r r 3 
B 2 1 



(3.15) 
(3.16) 



T ee = 3-[r-2Af-2(r-Af)cos 2 6l], (3.17) 

, [r-2M + 2Afcos 2 6»l. (3.18) 
87r r 3 sin 2 6> L J V ' 



This energy-momentum tensor is the source of the met- 
ric perturbation that will be calculated in the following 
subsections. 

It is easy to see from Eqs. (|3.12|) - (|3.18f) that the angu- 
lar dependence of the energy-momentum tensor is con- 
tained entirely in spherical-harmonic functions of degrees 
I = and 1 = 2; and because there is no dependence on 
(/), only functions with azimuthal index m — are in- 
volved. It can also be seen that the angular dependence 
of the energy-momentum tensor has an even parity. Our 
solution to the equations of black-hole perturbation the- 
ory will therefore have the following properties: (i) it 
will be axially symmetric; (ii) it will contain even-parity 
spherical-harmonic modes with (l,m) = {(0, 0), (2, 0)} 
only; and (iii) it will be stationary. The metric perturba- 
tion will contain a term of magnetic origin, and it will also 
contain a homogeneous term associated with the ambient 
Weyl curvature. 



2 (r - 2M)(r - AM) 



Q b = 2b 2 



(3.21) 
(3.22) 



where 2 b 2 := 2^B 2 /3. 

We now integrate the perturbation equations for the 
two relevant functions K (r) and h vv (r) — please refer to 
the listing of field equations in the Appendix. We first 
substitute Eq. into Eq. (|2~5f)|) and solve for K. The 

general solution is K(r) = —^b 2 r 2 +p + q/r, wherepand 
q are arbitrary constants. As discussed in Sees. II C and 
II E, we may exploit the remaining gauge freedom to set 
them equal to zero. We have, therefore, 



K 



2„2 



-Yr 



(3.23) 



The remaining field equations provide a number of equiv- 
alent differential equations for h vv . The general solution 
is h vv (r) = -i& 2 r(3r - 8M) + 25M/r. It involves an 
arbitrary constant 5M that can be interpreted as a shift 
in M, the black-hole mass parameter. To reflect the fact 
that we wish our perturbed black hole to have the same 
surface area as our original, unperturbed black hole (this 
was motivated back in Sec. Ill A), we set SM = 0. We 
will verify in Sec. Ill F that this condition does indeed 
lead to a preservation of the horizon area. We have, 
therefore, 



hyy 



1 



b 2 r(3r - 8M). 



(3.24) 



Sub stituting Eqs. lj3~2~3"|l and (|3~2l|) into Eqs. $2~Wfi - 
(|2~T2|) yields 



Pw = --B 2 r{3r-8M), 
Poo = ~~B 2 r\ 

= - -B 2 r* sin 2 6 



(3.25) 
(3.26) 
(3.27) 



for the I = sector of the metric perturbation. 



C. Integrating the perturbation equation: / = 

As discussed in Sec. IV D of Martel & Poisson [2lJ (see 
also the Appendix of this paper), the relevant projections 
of the energy-momentum tensor when I = are Q ab and 
Q b , which are defined in the Appendix. Using F 00 = 
1/V47T and the energy-momentum tensor of Eqs. H3.12JI - 
(|3.18|l . we obtain 



4& 2 , 



Q' 



, r - AM 



(3.19) 
(3.20) 



D. Integrating the perturbation equation: / = 2 



The relevant spherical-harmonic functions are Y 20 = 
\^/bJ^{2,cos 2 e - 1), Y 2Q = -^^/b/nsmecose, Y 20 = 



V20 _ 



/ 5/nsm 2 6, Y t 



0, and Y 



20 



— |^/5/7rsin 4 9. The required projections of the energy- 
momentum tensor are defined in Eqs. (4.17)-(4.20) of 



2 The constant b is not to be confused with the functions b (v) 
introduced in Sec. II C. 
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Martel & Poisson [21| (see also the Appendix of this pa- 
per); they are 



Q vv 


= -b 2 , 


(3.28) 




, 2 r-M 




Q vr 


= -b 2 , 

r 


(3.29) 


Q rr 


b2 (r-2M)(r-M) 


(3.30) 


Q" 


= -b 2 r, 


(3.31) 


Q r 


= -b 2 (r-2M), 


(3.32) 


cr 


= b 2 1 


(3.33) 




= -b 2 r{r-2M), 


(3.34) 



wh ere t he constant 3 b 2 has been reassigned to b 2 := 

We now integrate the field equations for the four func- 
tions h vv (r), j v (f), K(r), and G(r). The equation for K 
decouples, as was shown in Eq. (|2.36ll . and it involves the 
source term Q vv . Exploiting the remaining gauge free- 
dom to set all integration constants to zero, we take the 
solution to be 



K = -b 2 r 2 . 



(3.35) 



The remaining field equations form a set of coupled or- 
dinary differential equations for the remaining quantities 
h V v, jv, and G. These equations are easily decoupled 
by taking additional derivatives, and we easily obtain 
general solutions to the higher-order equations. These 
would-be solutions involve a number of integration con- 
stants that are not part of the true solution space; these 
are determined by substituting the would-be solutions 
into the original system of second-order equations, and 
making sure that the solutions stay valid. At the end of 
this process we obtain h vv = —c\M/r 2 — ic 2 (r — 2M) 2 + 
\b 2 M{r - 3M), j v = |ci(r + M)/r - c 2 r 2 (r - 2M) + 
\b 2 Mr 2 , andG= \c x jr - c 2 {r 2 - 2M 2 ) + \b 2 {r 2 + A/ 2 ), 
where c± and c 2 are the remaining constants of integra- 
tion. As we shall show below, the gauge freedom that is 
still at our disposal can be exploited to set c\ = 0. Setting 
also c 2 = i& 2 — for later convenience (thus discarding 
C2 in favor of the new constant e), our solutions are 



h vv = --b 2 (ir 2 - UMr + 18M 2 
6 



(3.36) 



Jv 

G 



+ e(r-2M) 2 , 

~-b 2 r 2 (r - 3M) + -er 2 (r - 2M), (3.37) 
6 3 



-b 2 M 2 + \e{r 2 - 2M 2 ) 
2 3 



(3.38) 



Substituting Eqs. l|335 | -<|3~38 ) into Eqs. ijO ty -pi ^ 



3 The constant b is not to be confused with the functions b m (v) 
introduced in Sec. II C. 



yields 



Pvv 



Pv9 



B 2 (3r 2 - UMr + 18M 2 )(3 cos 2 6 - 1) 

9 

+ £{r -2M) 2 (3 cos 2 6 - 1), (3.39) 
^B 2 r 2 {r -3M) sintfeostf 
- 2£r 2 {r- 2M)sin0cos0, (3.40) 
is 2 r 4 (3 cos 2 9 - 1) + B 2 M 2 r 2 sin 2 6 
+ £r 2 (r 2 - 2M 2 )sin 2 0, 



(3.41) 



v<t>4> 



-B 2 r A sin 2 0(3 cos 2 9 - 1) - B 2 M 2 r 2 sin 4 1 



£r 2 {r 2 -2M 2 )sin 4 



(3.42) 



for the I = 2 sector of the metric per turbation. We have 
introduced the constant £ := i^/5/7re; its interpretation 
as a tidal gravitational field will be examined below. 

We must now explain why it was admissible to set 
ci = in our solutions. We go back to Eqs. (j2.19jl (|2.22|l 
and consider the subclass of gauge transformations that 
leave K unchanged (in addition to h vr , h rr , and j r , which 
are all zero in the light-cone gauge). We see that when 
1 = 2, the subclass is characterized by a single function 
a(v), with the other functions related to it by b(v) = 3a 
and c(v) = — |ct. Taking a to be a constant produces 
c = 0, and we observe that under such a gauge transfor- 
mation, h vv changes by a term 6aAI/r 2 , j v changes by 
a term — 2a(r + M)/r, and G changes by a term —2a/r. 
We then see that selecting a — |ci produces a gauge 
transformation that effectively sets c\ to zero. There is 
therefore no loss of generality in making this assignment. 



E. Perturbed metric 

Combining Eqs. i|3~!^ - lj3~27|l from Sec. Ill C and 
Eqs. (|3~39V(|3~4*2l from Sec. Ill D gives us the metric of 
our perturbed black hole. Its nonvanishing components 
are 



9w = -f-^B 2 r(3r-8M) 
1 



-B 2 {3r 2 ~ 14Air + 18M 2 )(3cos 2 0- 1) 
9 



+ £(r-2M) 2 (3cos 2 0-l), 

9vr 1) 

2 

g vg = —B 2 r 2 (r — 3Af) sin cos 
3 

-2£r 2 (r-2M) sin cos 0, 

gee = r 2 - ^B 2 r A + ^B 2 r 4 (3 cos 2 - 1) 

+ B 2 M 2 r 2 sin 2 

+ £r 2 (r 2 - 2Af 2 )sin 2 0, 

= r 2 sin 2 0- -S 2 r 4 sin 2 
9 



(3.43) 
(3.44) 

(3.45) 
(3.46) 



9 



+ ^-BV sin 2 6»(3cos 2 ( 

-B 2 M 2 r 2 sin 4 
-£r 2 (r 2 -2M 2 )sin 4 l 



(3.47) 



and we observe that the perturbation is small when- 
ever r 2 B 2 -C 1 and r 2 £ <C 1, as was anticipated in 
Sec. Ill A. This is a three-parameter family of solutions 
to the Einstein-Maxwell equations, accurate through or- 
der (B 2 , £). The electromagnetic field is generated by the 
vector potential of Eq. I|3.3|) : it is accurate through or- 
der B. The parameters of the family are the black- hole 
mass M, the magnetic field strength B, and the tidal 
gravitational field £. 

The interpretation of £ as a tidal-gravity (Weyl- 
curvature) parameter comes from an examination of the 
asymptotic behavior of the metric when r> M (keeping 
r <C l/B, as was discussed in Sec. Ill A). In this regime 
Eqs. (E33h<E32l> reduce to 



-1 - -B 2 r 2 - -B 2 r 2 (3cos 2 9- 1) 
3 3 v 

f £r 2 (3cos 2 6>- 1), 



Qvr 



= l, 



(3.48) 
(3.49) 



-B 2 r 3 sin6»cos6» - 2£r 3 sin cos 0, (3.50) 
3 



r 2 --B 2 r 4 + -B 2 r 4 (3cos 2 0-l) 
9 9 y ' 



£r 4 sin 2 0, 
2 



(3.51) 



9U 



r 2 sin 2 



1 



9 



-B 2 r 4 sin 2 I 



g B 2 r 4 sin 2 0(3 cos 2 - 1) 
£r 4 sin 4 0. 



(3.52) 



The asymptotic metric no longer refers to the central 
black hole. It is the metric of a spacctimc that contains 
only a uniform magnetic field, expressed in an advanced 
coordinate system that is adapted to the incoming light 
cones of an observer situated at r = 0; the metric is 
limited to a domain r < a, where a is a length scale such 
that a 2 B 2 <C 1. The observer, of course, is fictitious, as 
r = is actually occupied by the black-hole singularity; 
nevertheless, the observer may be thought to exist in an 
unphysical extension of the asymptotic spacetime beyond 
its domain of validity, r ^> M. 

The metric of an arbitrary spacetime in light-cone co- 
ordinates was thoroughly investigated in our compan- 
ion paper [2{|. By comparing our Eqs. (|3.48|) - (|3.52|) to 
Eqs. (4.9)-(4.12) of the companion paper, we infer that 
the asymptotic spacetime is characterized by the follow- 
ing irreducible quantities: p := B 2 /(8ir) is the mass- 
energy density of the magnetic field as measured by the 
observer at r = 0, Su = S22 = — §£33 : = -B 2 /(127r) are 
the nonvanishing components of the tracefree part of the 
field's stress tensor, and T := B 2 /(8ir) is the trace of the 
stress tensor; these assignments are precisely what should 



be expected for a uniform magnetic field, 
ison reveals also that £u = £22 = — \8% 



The compar- 
:= £ are the 



--22 — — 2 033 

nonvanishing components of the spacetime's Weyl curva- 
ture tensor. (The irreducible quantities are all defined in 
our companion paper.) The comparison therefore gives 
us an operational meaning for the parameter £ : As was 
already anticipated, it is the Weyl curvature (the tidal 
gravitational field) of the asymptotic spacetime as mea- 
sured by an observer comoving with the black hole in the 
region M < r < 1/B. 



F. Perturbed event horizon 

The perturbed black-hole spacetime retains <j> a as a 
rotational Killing vector, and it retains t a = (1, 0, 0, 0) as 
a time-translation Killing vector. This vector is timelike 
outside the black hole, but it becomes null on the event 
horizon (which is therefore a Killing horizon). Setting 
g a /3t a t^ = g vv = and involving Eq. (|3.43(l informs us 
that the event horizon is now described by 



'■ n zon (6) : 2M( ! • -.\/'/;'M,r 



(3.53) 



It is interesting to note that JV>rizon(0) involves B 2 but 
not £. 

The horizon's intrinsic geometry is obtained by insert- 
ing Eq. 1 (535) 1 into the perturbed metric. It is described 
by the two-dimensional line clement 

dsLizon = 4M 2 [l + M 2 (B 2 + 2£) sin 2 0] d0 2 

+ AM 2 sin 2 8 [l - M 2 (B 2 + 28) sin 2 0] d<t> 2 . (3.54) 

The element of surface area on the horizon is 
4M 2 sin d0d(f>, and the integrated area is 



.4 



horizon 



16ttM 2 



(3.55) 



As was anticipated in Sec. Ill A, the perturbed black hole 
has the same surface area as the original Schwarzschild 
black hole; this reflects its quasi-static and reversible im- 
mersion within the magnetic field. 

The distortion of the event horizon can be measured 
by the Ricci scalar associated with the two-dimensional 
metric of Eq. (|3.54|) . This is 



R 



1 



2M 2 



1 + 2M 2 (B 2 + 2£)(3cos 2 0-1) 



(3.56) 



The distortion has a quadrupolar structure. The larger 
concentration of curvature at the poles reflects the 
greater strength of the magnetic field there; as was ob- 
served back in Sec. Ill B, the square of the magnetic field 
is given by |£?| 2 = B 2 cos 2 8. 

It is interesting to note that in accordance with the 
zcroth law of black-hole mechanics, the horizon's surface 
gravity displays no trace of this distortion. The surface 
gravity n is defined by the statement that on the horizon, 
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t a satisfies the generalized form of the geodesic equation: 
i^V 'pt a — nt a . A short calculation based on this equa- 
tion reveals that k — 1/(4M) plus terms of order B 4 , 
B 2 E, and E 2 . The surface gravity is uniform on the hori- 
zon, and it keeps its unperturbed, Schwarzschild value. 



G. Comparison with the Schwarzschild-Melvin 
solution 



There exists an exact solution to the Einstein-Maxwell 
equations that describes a nonrotating black hole im- 
mersed in Melvin's magnetic universe l35l|36ll37l. Known 
as the Schwarzschild-Melvin solution |22l l23l l24 | . it has 
a metric given by 



ds 2 = A 2 (-/ dt 2 + f- 1 df 2 + f 2 d9 2 ) + A~ 2 f 2 sin 2 ■ 



and a vector potential given by 



A° 



^BAcb a , 



(3.57) 
(3.58) 



where <f> a := (0,0,0,1) is the spacetime's rotational 
Killing vector. We have / := 1 — 2M/r as before, and we 
introduce the function 



A := 1 + -B 2 f 2 sin 2 
4 



(3.59) 



This is a two-parameter family of black- hole solutions; 
the first parameter is the black-hole mass M, and the 
second is the magnetic field strength B. 

The solution of Eqs. I|3.57|l - I|3.59|) is exact, and we 
wish to compare it with the perturbative solution of 
Eqs. (|3.43|) - (|3.47|l . We must first linearize the exact so- 
lution with respect to B 2 and transform the coordinates 
from the original system (i, f, 9, </>) to the light-cone sys- 
tem (v,r,9,<j>). The transformation from t to v is the 
same as for the Schwarzschild spacetime: 



v = t + f + 2Mln(f/2M - 1). 



(3.60) 



The transformation from f to r is designed to change the 
g vr component of the metric tensor from its current value 
A 2 ~ 1 + i_B 2 f 2 sin 2 9 to the new value of 1. It is given 
by 



1 + -B 2 f 2 sin 2 9 + 0(B 4 ) 
6 



(3.61) 



The angular coordinates (9, <f>) are not affected by the 
transformation. 

These manipulations bring the Schwarzschild-Melvin 
metric to the new form 

g vv = -/- -B 2 r(3r-8M)sin 2 + O(.B 4 ), (3.62) 



1 + 0(B 4 ), 

--B 2 r 3 sin 9 cos 9 + 0(B 4 ), 



1 



9ee 



T + -B' 2 r 4 sin 2 9 + 0(B 4 ), 
6 

r 2 sin 2 9 - -B 2 r 4 sin 4 9 + 0(B 4 ). 
6 



(3.63) 
(3.64) 
(3.65) 

(3.66) 



Comparison with Eqs. (|3.43|) (13.4711 reveals that the so- 
lutions are identical provided that we restrict the param- 
eter freedom of the perturbative solution. Indeed, to get 
a match we must set 



(3.67) 



The Weyl curvature of the Schwarzschild-Melvin solution 
is intimately related to its magnetic field. This feature is 
in fact inherited from Melvin's pure magnetic universe, as 
can be inferred from reading Sec. V B of our companion 
paper |25| . 

We conclude with the following statement: While the 
Schwarzschild-Melvin solution has the advantage of be- 
ing an exact solution to the Einstein-Maxwell equations, 
the perturbative solution of Eqs. H3.43|) - (|3.47|) has the 
advantage of possessing a larger number of parameters. 
The perturbative solution can therefore represent a wider 
class of physical situations. In particular, it provides 
the description of a magnetized black-hole spacetime in 
which the tidal gravity is not directly tied to the magnetic 
field. 
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PERTURBATION EQUATIONS IN THE 
LIGHT-CONE GAUGE 



In the even-parity sector the nonvanishing perturba- 
tion fields are h vv (v,r), j v (v,r), K(v,r), and G(v, r). 
According to Eqs. (4.13)-(4.16) of Martel & Poisson plf . 
they satisfy the field equations 
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or 

Q" 
Q' 



Q2 

Qy2 


P) 

--—K, 

r dr 




92 K 
dvdr 


2 9 TS 

r dv 


] 


dv 2 


r-M 


d 




dv 


ur + 


AM 


AC 



r dr 
K + 



2r 2 dr v 



— M d 



2r 3 



r dv 
M) 

— jv 



r 2 dv 



if/ 

2r 2 



A" 



4r 2 



9r 
/ 9 



A 



^2 



~3-> v 2r 2 4r 2 



r 

m)/ a 



dr 



dr 



K 



G. 



^-j v -± K -H—G--j v 
dr 2 v dr 2 dr r 2 v ' 



d 2 



dvdr 

d 2 



Jv 



Jv 



d 



K 



^9 nf d 



dvdr 
~2r 2 



K 

d 2 



dv 

d 2 



2 dv 
d 2 



G 



2d_ 

r dv 



-^-K+f-^K--^h v 
r dv dr z or z 



2 dr 



G 



29 h 
r or 



4 K+ 

dr 
A__d 

r 2 dr'' 



dr 



> Jv ! 



2(r-M) d 
r 2 dr 



K, 



dvdr 



G 



2rU 
dv 



2 9 2 

/ Q r 2 



r 



where A := Z(Z H- 1) = fi + 2 and ^ := (l-l)(l + 2) = A-2. 
According to Eqs. (4.17)-(4.20) of Martel & Poisson [2ll| . 
the source terms are 



8tt / T ab Y lm dn, 
16nr 2 



T aByVm dn _ 



1(1 + 1) 
8nr 2 J T AB n AB Y lm dfl 
32nr 4 



(1-1)1(1 + 1)(J + 2) 



where 



(f,r). The perturbation equations are not 



all independent; they are linked by the Bianchi identities 
d d M 

= ttQ vv + -^-Q vr + ~yQ vv 

ov Or r z 



-Q l 



Mf QVV 







_A_ r 1 

r " 2r 2 
d d 

TTW + TT 1 ^ "I 

cm or r 

+ -Q rr ~ AQ r ~-Q\ 

r Zr z r 

d d 2 

-Q r + -Q r + Q 



cm 



dr 



2r 2 



When I — the only nonvanishing perturbation fields 
are h vv and K, and the only relevant equations are those 
involving Q ab and Q b . When I = 1 the only nonvanish- 
ing perturbation fields are h vv , j v , and A, and the only 
relevant equations are those involving Q ab , Q a , and Q b . 

In the odd-parity sector the nonvanishing perturbation 
fields are h v (v,r) and h,2(v,r). According to Eq. (5.8) 
and (5.9) of Martel & Poisson [2l|, they satisfy the field 
equations 



d 2 



h 2 7rh v 



dr 2 v 2r 2 dr 

d 2 2d 
P = —7r-^-h„ + -Tr-h 



p = - 



dvdr v r dv v 2r 2 dv 

d 2 19. / d 2 

— — -h 2 H — ~£~hz — + 

Ovor r Ov l dr z 



H d fifd u uf 



2r 2 dr 
r_ZM_d_ } 

r 2 dr 



9 r - 4M 

Alt, o <l2, 



<9r 



r 



where A := 1(1 + 1) = n + 2 and \i := (Z-l)(Z + 2) = A-2. 
According to Eqs. (5.10) and (5.11) of Martel & Poisson 



1211. the source terms are 



P" 



16-Kr 2 
1(1+1) 



T aB x lrn ^ 
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167rr f ^ab j£i m ^ When I = 1 the only nonvanishing perturbation field is 

(I — 1)1(1 + + 2) J AB ' h v , and the only relevant equations are those involving 



The perturbation equations are not all independent; they 
are linked by the Bianchi identity 

o = ^-p v + ^P r + -P r - 4p. 

ov or r r z 



pa 
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